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This  report  provides  ell  the  necessary  Intonation  to  nake  practical  appli- 
cation of  acceptance  tests  based  on  a prior  distribution,  when  an  exponential 
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the  problam  and  a sunmary  of  the  other  voltaes.  Volume  II,  "Bisk  Criteria  and 
Their  Interpretation,"  defines  the  different  statistical  risks  which  can  be 
used  to  characterise  the  tests  and  provide  a guide  for  selecting  appropriate 
risks  la  various  practical  situations.  VoIixm  111,  "Implications  and 
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Determination  of  the  Prior  Distribution,"  provides  the  Mens  for  determining  the 
parameters  of  the  prior  distribution  from  existing  data,  and  discusses  the 
reason  for  using  an  inverted  gamma.  Volume  IV,  "Design  of  Testing  Plans," 
provides  instructions  for  establishing  a test  time  and  number  of  allowable 
failures  based  on  the  prior  distribution  and  the  selected  risks.  Volume  V, 
"Sensitivity  Analyses,”  shows  the  effects  on  the  test  parameters  caused  by 
changes  in  the  prior  parameters. 
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l.  nmcoucncN 

Consider  a sequence  of  electronic  oonponants  having  an  exponential 
failure  density  given  by  f (tl • ) - exp (-t/0 ) ; t > 0,  e > 0. 

Due  to  the  inherent  variability  in  the  production  process,  the  man  time 
to  failure  e,  which  characterizes  each  ocnponent,  will  differ  from 
ocnponent  to  ocnponent.  For  seme  production  processes,  the  variation  in 
6 may  be  so  small  as  to  be  practically  negligible  and  all  oonpenents  my 
be  assuwBd  to  have  a constant  0 . Other  processes  my  exhibit  a slowly 
drifting  tendency  such  that  oenpanents  produced  closer  together  in  time 
have  the  sane  e where  as  oonponants  produced  farther  apart  in  time  have 
different  values  of  0.  In  a very  general  setting,  the  sequenoe 
® , i ” 1*2, . . .can  be  considered  to  farm  a correlated  Hmp  with 

Poseible  ncn-staticnarity  and  seasonality  and  must  be  construed  as  one 
realization  of  the  underlying  stochastic  process  for  0.  One  such 
secpienoe  of  0's  is  shewn  in  Figure  1.  1.  In  this  report,  however,  we  will 
only  consider  the  case  of  independent  but  varying  0's. 

At  this  point,  it  is  necessary  to  distinguish  between  a ocnponent,  a 
system  and  a lot  in  the  context  of  reliability.  A ocnponent  is  a device 
vhich  is  not  repaired  upon  failure  such  as  a light  or  a vaocuun  tube. 
A system  consists  of  a large  lumber  of  oenponents  and  we  assure  that  upon 
failure,  can  be  repaired  to  be  as  good  as  new.  A lot  consists  of  a large 
nunber  of  oenponents.  Usually  oenponents  produced  within  a certain  tine 
span  are  groped  together  in  a lot. 

No#  we  discuss  two  types  of  models  for  describing  variation  in  0. 
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Figuee  l.lt  A WMlIftion  of  ttw  Gbcxolatad  stquanx 
tteOwMRta. 


For  the  cbm  of  a lot,  the  components  in  each  individual  lot  have 
different  e*,  i ■ 1,2,.. .drawn  independently  from 9(6).  The  lots 
are  aamanad  to  have  identical  9 (0) . This  case  is  illustrated  in 
Fig.  1.3. 

A generalisation  of  1.2  (a)  is  to  assure  different  9(6)  far  different 
lots  as  shewn  in  Fig.  1.4. 

Briefly,  population  heterogeneity  model  assures  that  each  population 
(lot  or  a system)  is  homogeneous  within  itself  with  a constant  e* 
and  e^'s  for  different  populations  can  be  described  by  a probability 
density  function  g(e) . Individual  heterogeneity  model  assumes  that 
each  population  has  a distribution  of  individual  6's  given  by  9(6) . 

In  either  case,  0 is  considered  to  be  a randan  variable  with  a 
distribution  g(  6)  called  tie  prior  distribution  or  the  oaipounding 
distribution. 


PRODUCTION  SKOUCNCt  OP  SYtTWi  OR  LOTS 


Figure  1.2;  A Asslizaticn  of  the  Indapenteit  of  0i,« 

for  System  or  Lots 


m this  report  m first  meins  ths  implications  of  ocraidering 
6 as  * randan  var labia  on  aana  aspects  of  rallability  analysis.  Speci- 
fically, tha  affact  of  tha  population  heteroganeity  on  aooaptanoa 

eanpling  and  tha  affact  of  individual  heterogeneity  on  failuza 

distribution,  aooaptanoa  sampling,  burn-in  and  systna  rallability  aza 
oonaidarad.  Tha  determination  of  tha  functioml  form  of  g(6)  aid  the 
reasons  behind  a specific  choloa  of  inverted  gene  prior  danaity  are 
diacuaaad  next.  Methods  for  determining  the  prior  paraneters  baaed 
upon  data  obtained  from  single  aeqale  truncated,  ceneozed,  with  re- 
pi  an—nt  and  without  replaoenent  aooaptanoa  aanpling  plana  for  system 
and  lota  are  given.  A detailed  analysis  of  aana  observed  failure  data 
is  presented  to  illustrate  the  process  of  obtaining  the  prior  density 


f. 


2.  SOKE  DffLICMIOWS  OT  PRIOR  IN  HELIABILnY 

In  this  Motion  we  consider  how  eons  of  the  reliability  analyses  are 
influenced  by  treating  6 as  a randan  variable  instead  of  a fixed  but  im- 
known  constant. 


2.1.  Population  Heterogeneity  Model 
Aoogptmoe  Sampling  and  Demonstration: 

In  reliability  acceptance  sanpling,  we  are  concerned  with  making 
aooept/re  ject  decisions  regarding  a sequence  of  lots  or  system  on  the 
basis  of  data  subject  to  random  fluctuations.  A definite  measure  of 
Iom  is  associated  with  each  of  the  two  decisions  when  they  are  inap- 
propriately taken.  Finally,  each  lot  is  sentenced  on  its  own  merit  without 
regard  to  decisions  taken  about  other  lots,  i.e.  the  sequence  ft1,  i ■ 1,2,... 
is  considered  to  be  an  independent  sequence.  In  reliability  demonstration, 
one  is  primarily  oonoemed  with  determining  whether  the  single  system  or  lot 
submitted  tor  testing  meets  the  specified  reliability  requirements  and  the 
concept  of  a sequence  of  lots  is  not  necessary. 

In  the  context  of  acceptance  sanpling  and  demonstration,  g(o)  may  be 
interpreted  as  a frequency  distribution  or  as  a degree  of  belief  distribution. 
The  two  interpretations  lead  to  philosophically  different  results.  A very 
brief  discussion  of  the  two  interpretations  is  given  below.  For  a detailed 
discussion,  the  reader  is  referred  to  Barnett  (1973) . 


» 
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Frequency  Interpretation: 

The  cl— leal  acceptance  sampling  plans  do  not  properly  consider  the 
relative  frequencies  with  which  different  values  of  8 occur.  This  dis- 
advantage can  be  removed  by  incorporating  g(6)  in  the  design  of  acceptance 
sampling  plans.  As  an  exenple  consider  the  design  of  single  sample  truncated 
plan  for  a system,  a system  is  put  on  a life  test  of  duration  T and  if  the 
observed  nunfcer  of  failures  r is  'ess  than -or  equal  to  the  acceptance 
nutter  r*  the  system  is  accepted;  otherwise,  it  is  rejected.  T and  r*  are 
determined  by  limiting  the  risks  of  wrong  decisions  to  a and  0 i.e. 
producer’s  risk,  P(R|e  ■ flQ)  i a and  consumer’s  risk  P(A|»  - 0^  $ 8,  where 
eQ  is  the  specified  value  of  «,  the  ndniaun  acceptable  value  of  6,  and 
A and  R denote  acceptance  and  rejection  respectively. 

This  approach  does  not  take  into  account  any  prior  information  avail- 
able regarding  0.  For  exanple,  if  it  is  known  a priori  that  P(0  - e^-ae)  * i 
than  all  systems  should  be  rejected.  However,  the  classical  (a, 8)  plan  will 
accept  100  0%  of  the  systems.  A solution  is  to  redefine  the  risk  criteria 
using  g(6).  As  an  exanple  consider: 

pr°du“r'*  »*•  *<■*>  v*«i«  «<•>  >»  - 

/P(A|0)  g(0)  d0  (1) 

Oonsixner's  Ride,  PfOstjA)  - - 


P (A) 


■ Winter  of  acoepted  systems  with  8^8, 
tlurber  of  accepted  systems  1 " b 


(2) 


■hHk 


Ths  risks  P(R)  and  0*  are  defined  in  the  long  run  average  sense  giving  due 
'•fcight  to  the  frequency  of  occurrence  of  8.  in  this  case  if  PO^)  = 1, 
then  PtesejA)  - 1 and  for  any  specification  of  (P(R),e*),  a test  plan  does 
not  exist.  Similarly  if  Pte^)  $ $*  then  POjjejA)  < 6*  and  all  the 
system  are  accepted  without  testing. 


Dcgre*  of  Belief  Interpretation: 

Consider  the  i system  under  acceptance  test.  The  system  has  0 ■ 0^ 
*iich  is  an  unknown  but  fixed  constant.  In  the  Bayesian  fonruLaticn, 
information  regarding  01  nay  be  quantified  by  the  subjective  prior 
g(0i) . If  the  only  available  information  is  that  Qi  is  independently 
dravn  from  g(0) , iters  g(0)  is  the  prior  distribution  based  on  a frequency 
interpretation,  then  g{0i)  = g(0).  However,  in  general,  gft1)  p g(6). 

In  fact  the  notion  of  g<0)  as  a frequency  distribution  is  dispensable 
and  one  may  use  a conjugate  prior  to  represent  g (0*) . 

If  the  system  is  put  on  a life  test  and  the  failure  tines  tj_,  tj...t 
are  observed,  then 


f(V  t2***tJel)“  “Tn  l Vei> 


» than  ragazdad  aa  a fraction  of  I1,  it  ia  kmn  as  tha  likelihood  fraction 

i J 

of  81  and  oontaina  all  lnfonetlon  regacding  • contained  In  tha  aanple 
i.e. 


then  flt^.y-  j L f(t1...thlii)  gd1)  dl1 
8 


1 .-(mi+l)  -<  ? VY,/8i 

gO^t,...!)- -li e1  e Jml 

1 71  r (n+x) 


(8) 


Equation  (8)  represents  the  density  of  an  inverted  game  distribution,  i.e. 
g(6i|t1...t^)  <v  Ga'  n+X). 


Die  design  of  a sapling  plan  can  be  acccnplished  by  confuting  appropriate 
posterior  risks  or  by  minimizing  a suitable  cost  function  using  the 
posterior  density. 


Differing**  in  the  Two  Approach* a: 

The  above  two  approaches  are  philosophically  different.  The  first 
anploys  g(0)  to  define  long  run  average  rides  and  is  meaningful  only 
ehile  dealing  with  a large  sequence  of  systems  or  lots  (eooeptanoe 
sapling) . The  second  considers  the  prior  as  a subjective  distribution 
and  the  approach  is  meaningful  in  a Bayesian  sense  whether  one  is 
dealing  with  a single  system  (demonstration)  or  a sequence  of  systems 
(eoosptanos  sapling). 


2.2  Individual  Heterogeneity  Model 
FcuZ&it  ViA&Ubution: 

Under  the  assumption  of  individual  heterogeneity,  a lot  consists 
of  components  with  different  e1,  i - 1,2, ...drawn  independently  from  g(e). 
If  a component  is  randomly  selected  from  this  lot,  its  life  tine  dis- 
tribution f (t)  is 

f(t)  - /f (t|e)  g(6)de  ( 9) 

e 

Assuming  an  exponential  conditional  failure  density  and  Ga'  (Y,X)  prior, 
we  have 


i.e.  t Pareto  (Y,X).  Hence  the  observed  lifetimes  tj^,  tj...  will  form 
independent  sstples  fran  Pareto  (Y,X)  distributions. 

Same  properties  of  f (t)  in  Equation  (1C)  are: 


E(t)  - r/(*-i)/  * > i 

(11) 

Var(t)- XY2/(X-i)2  (X-2);  X > 2 

(12) 

(13) 

h(t)  - X/h+t) 

(14) 

Note  that  Pareto  has  a decreasing  hazard  rate  whereas  exponential  has 
a constant  hazard  rata. 

Itus,  if  the  individual  heterogeneity  model  applies,  the  acceptance 
sanpling  plans  ahculd  be  based  upon  the  Pareto  distribution  rather  than 
the  exponential  distribution. 

8 u/in-in 

It  oonponents  fran  a heterogeneous  lot  are  put  on  a bum- in  test  of 
duration  T,  the  weak  individuals  (oonponents  with  small  6)  will  be 
eliminated  and  the  population  hazard  rate  will  decrease.  Specifically, 
from  Eq.  (14' , h(0)  * \/y  and  h(T)  ■ X/ (y+T)  . With  individual  hetero- 
geneity, bum-in  will  inprove  reliability.  For  the  case  of  population 
heterogeneity,  burn-in  has  no  desirable  effect. 

Let  us  suppose  that  the  oonponents  are  required  to  meet  a reliability 
requirement  R for  a mission  time  T^.  The  necessary  bum- in  time  may  be 
computed  as  follows: 

, R(TfTl)  v +T  > 

ROM.  |t  > T)  — - (J^or-)  - R (15) 

l R(T)  Y+T+t1 

(Y-HT, ) R^-y 

Henoe  T = (16) 

1 - Rr' A 

The  expected  fraction  of  parts  lost  during  bum- in  is 

1 - R(T)  - 1 - (-L-)X  . 
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Optinun  bum- in  time  can  be  determined  on  the  basis  of  cost  consider- 
ations. The  oost  nodal  will  include  the  oost  of  bum- in  (use  of  bum-in 
facility  and  the  oost  of  ocnponents  lost)  and  the  savings  due  to  in- 
creased population  reliability. 

St/Atem 

Consider  the  case  of  n ocnponents  in  series.  The  system  is  un- 
successful if  ary  one  of  the  ocnponents  fails.  Under  the  usual  exponential 
assumption  the  system  reliability  is: 

R_  (t)  - fl  R.  (t)  = n e"^®  - e~nt/®  (17) 

3 i-1  1 i-1 

However , for  the  case  of  population  heterogeneity,  with  an 
inverted  gamna  for  0,  the  system  reliability  is  given  by 

V«  - | - (yfet)  * <“> 

For  the  case  of  individual  heterogeneity,  the  system  reliability  is 

R At)  - n R.  (t)  - n (-4r)  - (rrr)1*  (19) 

8 i-1  1 i-1  ***  Y+t 

Zf  there  are  N different  types  of  ocnponents  and  there  are  n^  ocnponents 
of  the  i**1  type,  then 

N y,  X-n. 

R«(t)  - n (--ip)  1 1 (20) 

8 i-1  V* 
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ya  and  Xi  ar»  the  parameter*  of  the  distribution  of  01  far  the  ith 
type  of  ocegmant. 

For  a parallel  system,  system  failure  occurs  if  and  only  If  all 
consonants  fail,  filiations  corresponding  to  (17),  (18),  (19)  anl  (20) 
are  given  by  the  following  Equation*  (21),  (22),  (23)  and  (24)  reapectively. 

Rs(t)  " 1 " * (1  “ \<t))-  1 - (l-e“^*)n  (2i) 

%<t)  - 1 - | (l-e”t/®)n 

R.  (t)  - 1 - 11  - (-*-)*  )n,  (22) 

■ Y+t 

Similar  expressions  may  be  obtained  for  other  types  of  systems. 


3.  CETHtCPWriOW  or  PRIOR  DENSITY 


The  choioe  of  the  functional  form  and  paxanatar  valuaa  of  the  prior 
density  is  critical  since  a wrong  choice  of  prior  will  seriously  vitiate 
subsequent  analysis.  Information  regarding  the  frequency  distribution 
g(0)  any  be  available  from  the  following  sources. 

• Field  (operational)  failure  data 

• Failure  data  fron  previous  acceptance  tests. 

• Failure  data  cn  similar  ocnponents  or  syetsme 

• Subjective  evaluation  based  qpon  a knowledge  of  system  design, 
failure  data  from  design,  development,  prediction,  assessment  and 
d— mat  ration  phases  etc. 

Subjective  information  regarding  the  frequency  distribution  g(6)  may  be 
quantified  using  the  'betting  odds'  approech.  Alternately,  if  the  form  of 
the  prior  density  is  known,  the  parameter  values  nay  be  'guMsed*. 

The  problem  of  empirical  determination  of  prior  is  somewhat  involved 
sinae  the  random  variable  • cmnt  be  directly  observed.  The  failure 
data  from  which  the  prior  must  be  obtained  is  usually  aveiljfcla  in  the 
following  two  format  for  a system  the  date  (T^  *A),  i - 1,2,. ..n,  are 
. available  there  TA  and  *A  repraeant,  respectively,  the  duration  of  life  test 
and  the  observed  matter  of  failures  for  the  1th  wrmtm  end  n is  the  mater 
of  similar  qrstame  on  s lifo  tost,  for  s lot  (ni#  TA,  j^),  i - 1,2,. ..N 
are  available  where  N la  the  mate  of  lota,  nA  la  the  mate  of  ocapcnanta 
from  the  i lot  an  test  for  duration  TA  end  is  the  matter  of  l'xrr+ 
in  Ta.  In  each  cess,  individual  times  ip  failure  nay  be  available. 
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Thus,  ths  aval  labia  data  generally  consist*  of  observations  on  sans 
random  variable  Y,  which  may  be  the  nutter  of  failures,  sample  NOP, 
time  to  failure  etc.,  and  not  from  the  distribution  of  0.  Hence  the 
marginal  distribution  of  y, 

f(y)  * /"f(y|0)  g(e)de  (25) 

o 

must  be  used  to  estimate  g (0) . 

It  should  be  noted  that  g(0)  is  a oonpountong  distribution  in  the 
oosmonly  used  terminology  (Ord,  1972).  However,  in  the  current  context 
g(6)  is  not  only  the  oospounding  distribution  but  also  represents  our 
prior  knowledge  about  6.  Mathematically,  there  is  no  difference  and 
results  of  ooapound  distributions  apply.  However,  to  properly  reflect 
the  role  of  9 (6) , it  is  called  the  prior  distribution  of  6. 

Oppose  9^(8)  and  g2(0)  are  two  priors  for  0.  For  using  f (y)  to 
estiamts  9(0)  it  is  necessary  to  ensure  that  if  the  following  relation- 
ship holds, 

• m 

f(y)  - / f(y|e)  g1(8)d0  -/  f<y|e)  g (*><*  (20 

o o * 

than  91<0)"  92(®>» 

This  idantlf lability  of  prior  is  crucial.  In  this  study  it  is  ensured 
by  results  dus  to  Hsichar  (1960,  1961). 

QLvmn  tha  failure  data,  ths  problem  of  estimating  g(0)  is  very 
ocapllcated  if  ths  fiann  of  tha  prior  density  is  unspecified  and  the  data 
acmes  from  unplanned  experiments.  In  the  following  we  assume  that 
f(t(0)  is  negative  exponential,  9(0)  is  inverted  game  with  parameters  (y,x) 
i.e.  g(0)  ' Ga*  (Y,X)  and  ( y,X ) are  estimatad  from  planned  aaparinenta. 
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3.1.  Inverted  G— » Distribution 

The  inverted  gna  prior  density  is  given  by 

9(4)  - ^ e"a+l)  e‘Ve  y»A#0  > 0 (27) 

g(0)  is  chosen  to  be  of  this  farm  because  (also  see  Schafer  et.  al.  (1970)): 

(a)  Subsequent  data  analysis  shows  Ga'  (y,A)  to  be  an  adequate  prior. 

(b)  As  seen  from  Figures  2.1  and  2.2,  inverted  gamma  is  very  flexible  and  is 
capable  of  representing  many  practical  situations. 

(c)  Being  a aon jugate  density,  it  ispliee  mathematical  ease. 

Inverted  ganma  is  a skewed  distribution  and  all  moment  \p  to  [A] 

The  median  and  mode  always  exist.  Seme  of  its  properties  are: 

Mean,  up  - y/(A-1) 

Variance,  op2  - y2/{\-l)2  (A-2) 

Mode,  y/  (k+l) 

3.2.  Parameter  Estimation  from  Planned  Experiments 
Me  now  consider  the  estimation  of  the  prior  parameters  y and  A based 

on  failure  data  obtained  from  planned  experiments.  Specifically,  failure 
data  arising  from  truncated,  censored,  with  replacement  and  without 
replacement  single  sample  plans  for  systems  and  lots  are  considered.  In 
each  case,  the  problem  is  reduced  to  one  of  fitting  the  appropriate  marginal 
densities  to  the  abeerved  data.  The  marginal  densities  for  the  various 


exist. 

(28) 

(29) 

(30) 


cams  wi  developed  in  the  following  subsections.  A description  of  the 
step  by  step  epproech  to  fitting  is  given  in  Appendix  A. 

Note  thet  if  e large  nuaber  of  observations  are  available  on  a large 
fiafcer  of  systems,  g(8)  can  be  directly  obtained  from  the  observed  values 
of  tines  to  failure.  Further,  g(9)  is  a consistent  estimator  of  g(6)  and 
hence  can  be  used  in  lieu  of  g(6). 

3.2.1.  Data  Fran  Truncated  Plan  For  a System;  x Type  II  Censoring) 

In  this  caae  n identical  system  are  put  on  a life  test  of  duration 

r. 

T each.  The  random  variable  X denotes  the  rnxtber  of  failures  in  time  T. 
The  times  to  failure  far  each  system  are  noted  and  the  resulting  data 
is  of  the  following  fount 


e1  e2  eA  + 

hi  hi  hi  h& 

hi  hi  fci2  fcn2 

e • e e 

e • e e 

e e e a 

X X X X 

\ \ \ \ 


Thus,  for  system  i with  an  unknown  MUST  the 


tine  intervals  t^,  ti2/**>#tixi  and  *** 


failures 


occur  at 

mm 

failures  is  t^. 
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Since  we  have  assumed  that 


f(t|6)  - ! , t $ 0,  « > 0, 

e 


it  follow*  that 


f(x|e)  - 


.•T/e/™/.,x 


The  marginal  distribution  of  x is 

90 

f(x)  - / f(x|»)  g(e)dfl 
o 


- • rk 

* t*  (Tfy)/e 

£<x)  " T WxT^  ps-xvi""  d0 

O 0 


By  letting  Z « (Tty)/e,  we  get 


f U)  - JC L_  /V2  z^+x”1^dZ  (34) 

-X*  TCI)  CW-A)a+x  o 
The  integral  tsxm  is  r (X+x) , whence 

f(x)  " mix'  x- 0,1,2...  (35) 

1 K Tty  Tfy 

i.e.  x has  a negative  binomial  density  with  paramater*  y,\.  Therefore 
estimates  y,  £ of  y and  X can  be  obtained  by  fitting  a negative  binomial 
distribution  to  the  observations  x^  Xj,...,^.  A catfxiter  pcogrm 
for  obtaining  such  a fit  is  described  in  Ascendix  A. 


3.2.2.  Data  From  Censored  Plan  For  a 


In  this  case  n identical  systems  are  put  on  a life  test  and  testing 
is  continued  until  each  system  experiences  r failures.  The  resultant 
data  is  of  the  following  form: 


m mmm 


- ’ ■ 
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1 


e1 

a2  .. 

e1 

ease  ° see 

...  tn 

hi 

hi 

hi 

hi 

ha 

e 

h2 

e 

h.2 

e 

h2 

e 

a 

s 

e 

hr 

a 

hr 

e 

hr 

Jl  r 
6 

0*-  F 
n 

Thus,  for  systam  i with  an  uiknom  MOP  9i  the  v failures  occur  at 

tins  intervals  t^,  t^,...^  and  the  average  ties  between  failures  is 

t, . bet  the  randan  variable  T denote  the  total  tine  an  test  for  a system. 
1 r 

T takes  values  T.-  \ t, .,  i - 1,2, ...n.  Since  T is  a suomation  of  r 
1 j-1  13 

exponential  randan  variables,  it  has  the  following  gamma  density. 

f (T|t)  • . e*^®,  T > 0,  e > 0 (36) 

P0r)er 

Using  the  prior  g(0)  fron  Bq.  (27),  we  gat 

ten  .flL  ,-<r*,)/t  ao  (J7) 

rcr)  rex)  ° 


l 
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By  letting  Z ■ (Tvy)/e,  we  get 


-£i  -xi. 


r(r)  r(X)  (Tty )' 


/ zx+r_1  e'Zd2 


The  integrated  term  is  r (X+r)  and  hence 


Let  y ■ T/y,  then 


- LktEL_  ^l1.! 
r(X)r(r)  CT*y) 


f(y)  -L&iEL-y^1  (l+y)”(x+r);  y > 0,  (40) 

rcortr) 


i,e.  y has  an  inverted  beta  distribution  with  parameters  X and  r or 
y * Be'  (X,r).  Sinoe  T/r  ■ (y/r)y,  we  get  • * (y/r)  Be*  (X,r).  Thus 
a rescaled  inverted  Beta  distribution  can  be  fitted  to  the  observed 

Ai  A i A.. 

values  6 , * to  get  estimates  of  y and  X. 


3.2.3,  Data  Fran  Censored  Plan  for  a Lot  (With  and  Without  Replacement) : 
In  this  case  k components  from  each  of  the  n lots  are  put  an  a life 
test  and  testing  for  each  lot  is  terminated  as  soon  as  exactly  r $ k 
failures  are  observed.  If  T denotes  the  total  time  on  test  for  a lot  then 
it  follows  from  Epstein  (1960)  that,  both  with  and  without  replacement 


2T/e  * *2  . In  the 
r * 

T « t ♦ Orrjt^.  In  each 


T ■ ktr  and  in  the  non-replacement 
» - T/r.  If  y - 2T/e,  then 


*^•1  ’ •‘1r/*  ’ *>»' 


(41) 


Hence 


fCr|») 


r(r)e 


i-  T^1  SV*,  T > 


(42) 


J* 

\ 

I 

I 


whirfl  *•  **•  ■**  gwma  density  as  given  by  Bq.  (36) . It  follows  that 


9 % h/t)  Be1  (X,r)  distribution  or  (£)3  has  an  inverted  Beta  distribution 


with 


X and  r. 


3.2.4.  Data  From  Truncated  Plan  for  a lot  (With  teplaoeramt) : 

^ this  case  K items  from  each  of  the  n lots  are  plaoed  on  a lift* 
test  of  duration  t^«  let  , Xj,,,,  denote  the  nunber  of  observed 
failures  in  lots  1,  2,  ...,n  respectively.  It  is  that  in  this 

the  conditional  distribution  of  tits  nvnbar  of  failures  in  t is: 


fH<» 

Hence: 


e"***/® 


f<*  . (A*  (%_)x  m.,  Til  , 0, 

r(x)jd'Kt_V  'Kt^V 


(43) 


(44) 


O ' ’ ~o 

which  is  negative  binoedal  with 


y and  X. 
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3.2.5.  Data  From  Truncated  Plan  tor  a lot  (Without  Repinf— w*) . 


in  this  csss  k itmm  tram  sach  of  the  n lots  are  plaoad  on  a life  test 
of  duration  tQ.  Let  R « 1 - axpf-t^e)  denote  the  reliability.  Osi 
(1-10  rayreeants  the  probability  of  failvme.  if  r denotes  the  roster  of 
failures  ten 

f (r|R)  - (J)  (1-R)r  (R)n"r  (45) 


Prior  for  R: 

Given  6 % Ga*  Cr,l)  and  R » e“^®# 


g(R}  --I e"(r“to)*'e 


r<M  t 


w 

I -1 

Since  e“(Y“to,^e  ■ R ° and  e-  - t^inR,  we  get 


\ 


(t/t  ) X , , Y/t-1 

g(R)  — C-inl)*  1 r ° 

rex) 


(46) 


(47) 


-t  A' 

If  R - e 511  ten  *nR  - (1-R)  and  if  , - Y/tQ  ► 0,  ten 
ax  » r(an)/  r(a),  for  email  i. 

This  gives  t 

9(W"  TCTFkr  ROrl  •”1*  0 < * < 1 (48) 

The  approximate  marginal  density  of  r is  thin  given  by: 


f to  - / f Cr|R)  g(R)dR  - (")  HatiJ . L<r»*-z)  r(rfi)  (49) 

O T r(a)  rft)  r(nfo+r) 

which  is  the  beta-binanial  distribution. 


4.  FllTlNG  COMPLETE  AND  TWUNCfllED  NEGATIVE 
BINOMIAL  DISTRIBUTION 


Ihe  purpose  of  this  section  is  to  illustrate  the  estimation  of 
prior  parameters  and  their  varianoe-covarianoe  matrix  for  the  case 
in  Section  3.2.1.  In  this  case,  the  data  consists  of  field  failures  on 
different  types  of  systems.  For  each  system,  the  observed  failure 
frequencies  nJt,  corresponding  to  the  observed  number  of  failures  x in 
a fixed  time  T,  are  available.  For  this  situation  it  has  been  shown  in 
Section  3.2.1.  that  ifs 


„-<*♦!>  e~r/e  e,Yjl  , 0 


fltle)  ■ t > 0 

*tei«-*2!{S2iL  X-  o,i,2 


'l„\  _ (A+X-l)  l , T .*  . V . . . 

f W ~77  ( ')  (-1-)  t X = 0,1,2,...  (52) 

xt  (A-l) ! Tty  Tty 


^rtheraare,  as  mentioned  earlier,  because  of  the  identifiahility  of 
the  prior,  if  X has  a negative  binomial  distribution  with  parameters  A 
and  > given  by  Bq.  (52)  and  the  conditional  failure  density  is 
exponential  then  g(0)  is  an  inverted  gamma  distribution  with  para- 
mstsrs  y and  A as  given  by  Bq.  (27) . 


Dw  parameters  (y,A)  can  be  estimated  by  fitting  a NBD(y,X)  to 
the  observed  data  (x,r^) . In  this  case,  when  the  failure  time  is 
exponential,  the  adequacy  of  the  negative  bincmial  fit,  which  is  a 
compound  Poisson  distribution,  ensures  the  adequacy  of  the  inverted 
gemma  prior  (See  Tsicher,  1960,  p.  71) . 

Sometimes  frequency  of  zeros  is  not  recorded.  In  that  situation, 
X follows  a truncated  negative  binomial  distribution  given  by: 


f(x>  - — x » 1,2, 
x*  (x-i) : {i  - ^ 


(53) 


Tty 


We  now  oonsider  parameter  estimation  for  the  oanplete  and  the  truncated 


negative  binomial  distributions. 

4.1.  Oonplete  Negative  Binomial  Distribution 

In  this  case  X has  a NBD(y,x)  distribution  given  by 


f (X)  - {■A+y-1!!-  (-Z- )X  x = 0,1,2,. 

x!  (X-l)  J ^Y  TfY 


with  mean  and  variance  given  by 


E(x)  - v - ^ , 

and 

VrrOc)  - 

y 


(54) 

(55) 
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There  is 


(p,X>. 


i advantage  in  estimation  if  we  write  f (x)  in  of 
wa  gat 


fcx)  - tesflg  (jl.\x  t±y, 

x!&-l)J  lyf*/ 


X - 0,1,2, ... (56) 


We  shall  consider  estimation  for  both  pairs  of  parameters  (y,x)  and 
(u»X). 

4.1.1.  Moment  Estimates: 

Equating  the  sanple  mean  x and  simple  variance  S2  to  the  population 
mean  and  variance  given  by  Bgs.  (54)  and  (55) , we  get 

y'  • x 


y'  - 


and 


XT 

S2-x 


(57) 


S2 

A'  - * , 

S -x 

kdiere  T is  the  test  time,  and  u' , y',  X'  represent  the  calculated 
values  of  u,  Y and  X,  respectively,  x as  an  estimate  of  u is  fully 
efficient.  It  should  be  noted  that  if  S2  < x,  the  negative  binomial 
may  not  be  appropriate.  (Also  see  discussion  on  page  60.)  It  is 
sham  by  Anscomb  (1950)  that  the  efficiency  of  the  moment  solution  is 
at  least  904  under  the  following  conditions: 

(a)  u is  saall,  £ or  y>6T 

(b)  is  large,  X>13 

(c)  u is  meditai,  (X-m)  (X*2)/u  > 15  or  (T*Y)  (X+2)/T  > 15 
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4.1.2.  Maximum  LlJoellhood  Estimates* 

Estimation  of  Y.X.* 

If  x takas  values  0,1,2. ..K  with  the  oortaspondong  observed  frequencies 
nQ,  n1#  n2>..nK  the  log  lUcalihood  can  be  written  as 


or 


Iat 


and 


Then 


or 


In  L - f n .In  f (x) 
x-0* 


(58) 


In  L - l n { (ln(X+x-l)!  - ln(x)!  -ln(A-l)! 
x-0  x 


♦ Xln(-JL)  + x ln(JL)} 

Tfy  Tty 


K 

l \ 

x-0  * 


nx  - l x n 
xH) 


ilnL  m nftlhrx) 
3y  (T*y)y 


22&-  11^  i - 

ax  »■!  j-l  (A+j-l) 


1 ■}  - n ln(l-f^) 


K 


" i.<". A.  * X v ■ n 111  (1^ 


ax  j-i  (x+j-i) 
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K 


(59) 


(60) 


(61) 


(62) 


(63) 


The  maximum  likelihood  estimate*  are  obtained  lay  solving 

. ainL 

»Y 


0, 


and 


31nL 

31 


- 0 


The  ten  simultaneous  equations  are  oomplicated  and  the  solutions  y and  x 
can  not  be  obtained  directly.  A general  method  is  to  assure  a trial 
solution  (y',X»),  which  may  be  taken  as  the  moment  estimates,  and  derive 
linear  equations  for  email  additive  solutions.  Such  a netted  of  scoring 
is  now  described  (Rao  1965) . 

Expanding  31nl/3y,  31nl*/8A  about  Cy  * # A 1 ) and  retaining  only  the 
first  power  of  6y  - y-y',**!  - X-X'  i~r\.  to 

Mi  *{Mi  + &y£}*L  + yX£lnL}  | - o 

3y  3y  By7  3y3X  I (Y-y»,X«X') 

+ + v^lsk)  i . o 

3*  3X  ax2  36 ax  I (Y-Y',X-X') 

The  incremental  corrections  are  obtained  as: 


(64) 


"«Y 

"32lhL 

32lnL~ 

- . 

3Y2 

3y3X 

6X 

32lnL 

32lnL 

" 3ydx 

3X2*^ 

L 

31nL 

3y 


I 8lnL  | 

L 3X  -J 


(65) 


| (Y-Y’») 


,X-X') 

The  parameter  values  after  the  first  iteration  are  Y"  ■ y'  + 4y 

and  X"  • A»  + &x.  H*  iterations  are  continued  until  oonvergence  occurs. 


second  partial  derivatives  In  Eg.  (65)  are  given  by 


2 — 2 
3 InL  m nxY  - nAT(Tf2v) 

y2  Cm)2  * 


3y* 


(66) 


3 InL  m nT 
3y9X  yCPty) 


(67) 


and  3lnL 


- lK  I !<-i  ?„) 

X-1  j-1  (X+j-1)2  j-1  (Xtj-l)2 


(68) 


Vca'ianct-CovcuPianoe  matrix  fop  y.^ 

The  varlanoe-oovarianoe  matrix  for  Y,^  is  given  by 


r Var(?) 


oov(y^> 


Oov(y,x)  var(x)  j 


_pa2lnL  E»2lnL  _~1 

3y2  ayax 


.A* 

r 

ay»x 


.32lnL 


ax 


(69) 


where,  from  Eqs.  (66) , (67)  arxJ  (68) 


Ei^.-nxX 

a7  y ttfY) 

gV.  ■« 

3y3X  y^y) 


(70) 


(71) 
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Estimation  oj 


Following  Filter  (1941,  1953) , Haldane  (1941) , and 


(1950), 


u and  x may  ba  obtained  as  follows.  The  log  likelihood  can  be  written 


Sira 


md 


aov(y.X)  » 0, 

as  MM 


(82) 

(83) 


var(x)  » - ■ . ■ ■ 

n(-Z-)2  1 ♦ 2 l 

j-2 


TTT(X)>i- 

(j-r 


(84) 


The  advantage  of  considering  (m,X)  m paraneters  radar  dan  (y.x)  is 
that  v and  X an  unoarralatad. 


jA.  Negative  Binoadal 

*•*"  nroa  an  not  recorded,  da  resulting  distribution  is  a tnneatad 
negative  binaeial  given  by 

*W  £i2EJ2i (-^-)X  (-I-)x,  x - 1,2....  (W) 

x!  (A-l) ! {1  - (JL)X)  **  T*r 
T*t 


The  mem  and  variance  are  given  by  [Johnson  & jfo t*  (1969)] 


L(x)  - (l  - (-*-)*>■* 


Tn 


\to(x)  - {]  . [1  - (J2,  ( 

y Tn  L frr  wt  J 


(86) 


(87) 


Equating  the  earpla  and  exported  values  of  mean  and  varianoe  gives 

«x>  simultaneous  equations  far  y and  X from  Wiich  mcment  estimates 

\ 

of  these  parameters  can  he  obtained.  The  equations  do  not  have 
explicit  solutions  and  Brass  (1958)  proposed  using  the  observed 
proportion  of  unit  values  along  with  x and  s2  to  obtain  y and  X. 

If  denotes  the  frequency  of  1,  than 

*T(n-n,) 

J (88) 


ns  - xtn-fij,) 


and 


X'  » 


nr'x  - n.  (Tn ') 


nr 


(89) 


T is  the  test  time,  y’  and  A' 

Y end  A respectively.  Brass  (1958)  also 
these  estimates  is  about  90«  in  most 
efficiency  of  thsse  estimates  is  given  in 


the  calculated  values  of 
that  the  efficiency  of 
A detailed  discussion  of  the 
4.2.2  and  4.2.3. 
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4.2.1  Maximum  Likelihood  Estlratss 


The  log  11  toll  hood  is  given  by 


InL  - l n {ln(x+x-l)J  - m(x)!  - ln(*-l)!  + xlnM-) 
*■1*  Tfy 


♦ xln(-3L)  - ln[l  - M-)  x] ) 

Tfy  Tfy 

on  differentiating  the  log  likelihood  function  we  get 


(90) 


ainL 


31nL 

TTJT 


nXT 


*T+y) 


nln  (=dr)  K 


nx 

7t+7T 


y — 7 ♦ L (7imT  L nx' 

l-^)  X j-1  U+3  x-j 


(91) 


(92) 


The  maxima  likelihood  equation*  obtained  by  setting  3lnl/ay  - o and 
31nL/3X  ■ 0 can  be  iteratively  solved  using  the  soaring  method  given 
earlier.  Penan  Lai  estimates  due  to  Brass  (1958)  nay  be  ueed  as  initial 


VurioriM-Covariano*  Matrix: 

An  estimate  of  the  variance-covariance  matrix  based  an  the  sarv-Je 
values  is  given  by 


Var(r)  cov(y,a) 


'_32lnL  _32lnL'1"1 

E 5—  E 

3y  3y  3 A 


.2 


Co v ( y » i ) var(A) 


3 InL  _3  inL 


1 e— — 
I £*3y3A 


where 


3 2 InL 


nAT  [(T+2y)-(T+2y+AT)  (^)  X]  nx 


y2  (T+y)  * }< 


(T+y) 


a 2 InL  nT  {1“<tIy)  ti-Alnt^j^-)  ] ) 


3yTX  y (t+y ) ( 


3 2 InL  _ n(T+Y)  {ln(f+Y)}  r r 1 r 

T v — I — TT 2.  I y l nx} 

3 A2  a-(-l_)A}2  j-1  (A+j-1)2  X-j 


A program  to  compute  Y and  their  varianoe  covariance  matrix  is  given 
in  Appendix  A. 


. . 


. .ata.jcirjriBVi  i [ 


4.2.2.  Efficiency  of  Estimation  Methods 

in  tills  ssctlon  w discuss  ths  efficiencies  of  various 
estimation  procedures  usad  for  estimating  the  paranetars  y and  X 
of  tiie  truncated  negative  binomial  distributions.  The  three 
methods  dlscusssd  in  ssctlon  4.2  are 

a)  Ths  method  of  ncomnts. 

b)  Ths  method  of  modified  moments  due  to  Erase,  and 

c)  The  method  of  maxima  likelihood. 

Ths  natural  question  that  arises  at  this  point  is  which  method 
is  ths  most  desirable  one  to  use.  m order  to  ampere  these  methods 
we  enploy  a widely  usad  criterion,  ths  Asynptotic  Relative  Efficirey 
(A.R.E)  given  by 


A.R.E.  ■ 


RejWPtotic  generalised  variance  x Information  determinant 
Ths  Aeynptotic  generalised  variance  is  given  by 
Var  (XJ  Oov  (X#,  yj 
COv  (A^,  y^)  var 


(97) 


lV“  «V  7.)|  - |_  - * | (99) 

ihave  X#  eed  y^  represent  the  estimators  of  X and  y respectively 
obtained  by  the  method  under  conparison. 

Ths  information  determinant  is  given  by 


^ ^ 1 l,j  " 1,2 


(99) 
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mm 


where  6^  - X,  Oj  ■ Y and  L Is  the  likelihood  fvncticn 
far  the  asynptotic  generalised  variance  for  various  a 
obtained  below. 


pressicra 

are 


a)  Method  of  Moeenta 

The  asynptotic  generalized  variance  for  the  manent  estimator 8 
X and  y la  given  by 

|Var  (X^,  ym)  | - lAj'  A2*  ^ Aj  A3I 
where  Aj,  Aj»  A3  are  as  defined  below. 

A^  is  the  variance  covariance  matrix  of  the  first  two  ample 

*3*  " "2  V) 

>V  - m2’2  / 

Aj  and  A3  are  the  Jacobdans  for  the  transfonnations 
V ')  (X,  p),  and 

(X,  p)  - (X,  y) 
respectively,  and  p - T/y 

For  the  truncated  negative  binomial  distribution. 


and  is  given  by 
*1 


/ V - V 
W " V **l' 


4.1) 


V 

‘*2’ 


i-{22» 

T 2 m 

X(X+1)  + Xffi 

ir«an'* 


(100) 

(101) 

(102) 

(103) 

(104) 

(105) 
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Define  Aj  ■ B_1.  Thm  Biia  2x  2 matrix  with 
,1,  »(I)  (2£lf\  log  <2SX> 


1(1.1)  - 


-X 


— — r 

<1  - (2£>  > 


(106) 


v2#T 


B(l,2) 


(2)  (2^) 


-x-i 


i-(2£) 


=X 


a - <£x>  > 


(107) 


l-  (2^ 


-X 


(2X+i)  (2)  + (2)  {x  (x+i)  (2)  + x<2)}(2Jl)  log  (2^1) 

B(2,l)  ■ " I. ™t  n 1 i ■■■  (108) 


Tfv-"1  2 

(l  - (^2)  } 


and 


2X  (X+l>  <2)  + X XJ  {(X+l)  + * 1 

b<2,2) ry 


<1  - <EI>  ) 


1-  flfft 

Tha  Jacobian  matrix  Aj  is  given  by 


ao9) 


1 

0 


(110) 


where  T is  the  test  time. 
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The  asymptotic  generalised  variance  can  be  obtained  by  plugging  in 
the  values  of  A^,  Aj,  Aj  into  Equations  (100)  and  then  A.R.E  is 
obtained  from  Equations  (97)  and  (99) . 

b)  Method  of  Modified  Moments  due  to  Brass 

Asymptotic  generalised  variance  far  Brass's  modified  moment 
estimators  and  i»  given  by 

|Var  ( X Ydt>)  | - Ia'5  A'4  AjI  (111) 

A4  is  the  variance-covariance  matrix  of  the  first  two  sanple 
moments  and  the  first  sample  frequency  and  is  given  by 

2 

1*2  - V V - V V 

A4-lu3*  -ux'  I i2*  u4'  -u2'2 

Ul-UjMPj^  (l-y2')Pl 

where  p^,  the  probability  of  one  failure  in  time  T,  is 

pl“ ^y)X  ^ (113) 

(i-  (^)  > 

Aj  is  the  Jacobian  for  the  transformation 

( Uj',  m2'  , Px)  (Y,X)  (114) 


a - Ul’)  \ 

(1  “ U2')  P 1 J (112) 

pl  ( 1 - P 
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the  elsoants  of  the  matrix  Aj  are  given  by 


*5(1,1)  - -T(.1"pf) __  . ui,T(1~Pi>  + a-p,)> 

v2  " Y2  " V (1"Pi>  <v2'  - Ux'  - u1'(l-pl))  2 


“ W,1  T (1-p.) 

Aj(2,l)  - i --1 

{w2‘  - u^2-  (1-Pj^) >- 


**«.»  - ^-^5 “L 2 T *y 

V2  ~ V * V d-Px)  (M2'  ~ Mj/2  - Wl*  (1-j^)} 


*5(1,2)  ■ (Ajd,!)  • + Y}/r 


*5(2,2)  ■ *5(2,1)  • uy  - Pl)/r 


*5(3,2)  » {*5(3,!)  • (y^  - Pl)  - (TfY)}/T 


c)  Method  of  Maximum  Likelihood 

Hie  Mynptotlc  generalized  variance  for  the  maxinun  likelihood 
•rtinetion  procedure  turns  out  to  be  V (Information  determinant) . 
nmefore,  it  is  clear  that  the  A.R.E.  of  MLE  is  alleys  one.  in  fact 
A.R.E.  Is  the  aspoptotic  efficiency  of  the  method  under  consideration 
OQB|>ared  to  the  method  of  likelihood. 


4.2.3.  A.R.E.  Carputaticna  and  Carpariaons  of  Estimation  Methods 

Um  asymptotic  relative  efficiencies  for  the  moment  estimators  far 
X ■ 0.1  to  5.0  and  p - 1.0  to  10.0  were  ocnputed  using  the  expressions 
given  in  Section  4.2.2.  and  are  given  in  Table  1(a)  • F*cm  this  Table 
we  note  that  far  constant  p,  A.R.E.  improves  as  X gets  larger.  Far 
ocnstant  X,  A.R.E.  deteriorates  with  increasing  p.  Thus,  A.R.E.'s 

are  generally  better  far  large  X and  snail  p values.  For  the  type 

\ 

of  data  that  we  have  analysed  in  Section  5,  we  find  that  X values 
are  nail  and  p mines  are  large.  The  A.R.E.'s  far  these  cases  are 
about  55%,  a rather  poor  value. 

A.R.E.  's  were  also  ocnputad  for  the  sane  range  of  values  of  X 
and  p for  the  method  of  Brass  and  are  given  in  Table  1 (b) . 

Me  see  that  the  A.R.E.'s  follow  a pattern  similar  to  that  of  the  method 
of  moments.  The  values,  however,  are  consistently  larger  for  the 
modified  moments  method. 

It  riuuld  be  noted  that  Braaa's  modified  moments  method,  vhich  is 
tw— a on  the  first  two  saqple  moments  and  the  first  sample  frequency 
is  computationally  such  simpler  than  the  method  of  momenta.  Also, 
a a seen  above,  A.R.E.'s  increase  considerably  by  using  the  extra  inform- 
ation extracted  from  the  flrst-saiple  frequency.  Even  though  the 
increase  in  A.R.E.  is  significant,  the  values  are  still  around  75%  for 
the  range  of  parameter  of  interest  in  Section  5. 
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lhs  minimum  chi-square  methods  described  by  Katti  and  Gurland 
(1962)  for  the  negative  bincmlal  distributions  can  be  easily  extented 
to  the  case  of  truncated  negative  binomial  and  they  will  certainly 
have  higher  A.R.E. 'a  than  the  moment  estimators.  These  estimators 
also  will  have  a simple  foam,  but  in  obtaining  these  estimators,  it 
is  necessary  to  solve  non-linear  equations.  Gurland  (1963) , and 
Hina  and  Gurland  (19(8)  have  described  the  generalized  minimum  chi- 
aquare  techniques  for  sane  generalized  Poisson  distributions.  These 
techniques  are  baaed  on  functions  of  ample  frequencies  and  ample 
moments.  These  procedures  lead  to  weighted  least  square  techniques 
and  the  estimates  can  be  obtained  as  solutions  of  these  weighted  least- 
squeree.  These  methods  also  can  be  easily  modified  for  the  truncated 
case. 

At  this  point,  it  should  be  pointed  out  that  the  generalized 
minima  chi-square  methods  will  have  higher  A.R.E.  when  mare  information 
is  extracted  from  the  data.  By  using  three  or  four  sample  moments 
and  two  or  three  sanple  frequencies  we  are  almost  guaranteed  (for  most 
of  the  two  parameter  families)  to  get  100%  efficiency.  However, 
these  sanple  functions  haws  large  aaynptotic  variances  when  higher 
momnts  and  frequencies  are  involved  and  henoe,  these  functions  change 
considerably  from  ample  to  sanple. 
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11m  method  of  maximal  likelihood  has  ease  good  asymptotic 
properties  and  this  method  works  fairly  well  in  oases  of  anall 
■■piss.  Also,  the  M£'s  always  have  a better  A.R.E  than  any 
other  generalized  minimal  chi-square  method.  Considering  all 
these  points,  w»  suggest  that  it  is  better  to  use  the  MLE  whenever 
it  is  possible  to  do  so.  Otherwise,  seme  generalized  minima  chi- 
square  method  ha  sad  on  three  ample  moments  and  may  be  two  or  three 
ample  frequencies  can  be  used. 

In  tiie  problem  that  we  are  considering,  it  is  not  difficult 
to  get  IC£'s  as  described  in  Section  4.2.1.  and  illustrated  in 
Section  5.  This  ensures  high  efficiency  of  the  estimators. 
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5.  PITTING  OF  PRIOR  DISTRIBUTICKS  FCR 


ekmt  a q*  sets 

W*  now  obtain  eetiastes  of  the  paraneters  and  their  standard  errors  for 
the  eight  sets  of  data  reported  lay  Schafer  (1970) . The  data  aonsists  of 
field  failures  on  eight  different  types  of  systems.  For  each  system, 
tl**  rbawsd  frequencies  nJ(/  corresponding  to  the  observed  nirter  of 
failures  x in  a fixed  Uirc  T,  are  available,  lhe  eight  data  sets  axe 
raproduoed  in  Tables  1.1  to  1.8. 

TO  further  explain  the  data,  consider  the  observed  frequencies  nx 
corresponding  to  the  observed  mnber  of  failures  x far  a search  indicator 
IP“128A  (Table  1.6).  fo  this  case,  a total  of  + 55  search  indicators 
woro  placed  cn  a life  test  of  duration  4320  hours  and  the  observed  mnber 
of  failures  for  each  search  indicator  were  recorded.  For  example,  21 
indicators  failed  once,  8 failed  twiae  etc.  The  frequency  nQ,  corresponding 
to  x ■ 0,  i.e.  for  the  indicators  that  did  not  fail,  was  not  recorded. 

Since  cfcwerviitions  corresponding  to  x ■ 0 are  missing,  a truncate*’ 
negative  binomial  distribution  aho  ld  be  fitted  to  the  observed  data  sets 
(x,!^).  Schafer  (1970)  however  fits  a oanplete  negative  binomial.  Table  2 
gives  the  predicted  frequencies  obtained  by  fitting  a truncated  negative  bino- 
mial distributions  to  the  date  in  Table  1.6.  For  ocapariscn  purposes  the 
predicted  frequencies  for  a oanplete  negative  binomial  distribution  are 
also  given  in  Table  2.  It  should  be  noted  that  the  aonplete  negative 
binomial  distribution  was  fitted  by  assuming  the  frequency  nQ*  0.  Table  3 
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• Schafer  (1970) 


TOTAL  NOSER  OF  FAILURES  74 


TABLE  1.4* 

OH  VIMO  JWggJFIBR  AM-472  (Data  Sat  4) 


HELD 


iulltJ. 


I 


TABLE  1.7* 

Fnqp  FMLUHE  DMA  ON  JWLIFIER  PANEL  PN-8  (Data  Set  7) 


XABE£  2 


OBSERVED  AtO  PBEMCIH)  FREQUENCIES  FOR  SEARCH 
INDICATOR  IP-128R 

Observed  muter  Observed  Predicted  Frequency  nx 

Truncated  Complete 

of  failures,  x frequency,  Negative  Binomial  Negative  Binomial 

______________  . 




0 

6.4 

1 

' 21 

19.4 

9.4 

2 

8 

11.0 

9.6 

3 

8 

7.1 

8.2 

4 

2 

' 4.8 

6^5 

5 

5 

3.4 

4.8 

4 

4 

2.4 

3.4 

7 

3 

1.8 

2.3 

8 

2 

1.3 

1.6 

9-19 

2 

3.8 

2.8 

55 

55 

55 

! 


58 


- w : ► . i :■  ?&r 


Ttotal: 


TORTf  3 


K2  qmPMESg  or  FIT  TEST  FOR  TOE  EKE*  OF 


Number  of  Observed  Ejected  A A a,2,a 

Failures  Frequancy  (n^)  Frequency  (n^)  Oye^)  v’W  /y 


Since  x2  flb|8lwl  - 24.2  and  X2  4 0#05  - 9.5,  aonplete  negative  binomial 


does  not  fit. 


Truncated 

Nvntoer  of 
Failures 


ive  Binomial  (role) 


Observed  Expected  A 

Frequancy  (n^)  Frequency  (y)  (nx-nx) 


<V$c>H 


Mill  55  55 

Sins  x2  , - l-5<  X2  2f0.05  - ‘.0.  «•  txvno.t-4 

Hfr"***1  is  a good  fit. 


contains  the  x2  goodness  of  fit  tests  far  the  fitted  distributions.  It 
is  clear  that  the  oonplete  negative  binomial  does  not  fit  whereas  the 
truncated  negative  binomial  distribution  indicates  a good  fit  to  the 
observed  data.  Moreover,  as  indicated  earlier,  the  oonplete  negative 
binomial  is  not  the  correct  distribution  to  fit  in  this  case. 

In  order  to  see  the  behavior  of  the  likelihood  functions  for  the 
data  being  analysed,  three  dimensional  plots  of  the  normalized  likelihood 
functions  (normalized  by  its  maximum  value)  for  data  sets  in  Tables  1.2, 
1.5,  1.6  and  1.8  are  shown  in  Figures  3.1  to  3.4,  respectively.  From 
these  plots  it  can  be  s^en  that  the  likelihood  functions  are  unimodal  in 
the  range  X>0,  y>0  and  thus  the  method  of  maxima  likelihood  employed 
here  will  in  fact  give  the  global  Hsudman  value  of  the  likelihood 
function  and  not  seme  other  local  maxima.  Using  the  method  of  maximum 
likelihood,  the  estimated  parameter  values,  asymptotic  standard  errors,  the 
observed  and  predicated  frequencies  and  the  x2  goodness  of  fit  tests 
for  the  eicfit  sets  of  data  of  Tables  1.1  to  1.8  are  given  in  Tables 
4.1  to  4.8.  In  each  case  a truncated  negative  binomial  distribution 
was  fitted  and  is  seen  to  provide  a satisfactory  fit  to  the  observed 
data.  The  results  are  simmarized  in  Table  5.1.  The  MI£  of  X for  Data 
Set  5 was  found  to  be  negative  and  we  have  taken  a value  0.00001  for  X 
in  this  case  as  explained  below.  Also,  from  Table  5.1,  it  is  inter- 
esting to  note  that  X far  all  eight  systems  lies  between  0 and  1.  In 

A 

view  of  the  standard  error  of  X,  X is  close  to  zero  which  large 

'variability  associated  with  the  mean  time  to  failure. 
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Thm  field  failure  data  an  Synch  Power  Supply  at»315  of  ^ir  1.5 
9ivea  a negative  value  far  the  MI£  of  the  permwter  A,  doaa  not 

belong  to  ita  acceptable  range,  m such  a caae,  a rnaenrrtili 
is  to  take  a proper  limiting  form  of  a given  distribution  giv 
a » 0.  However,  die  value  A - 0 is  also  not  aamptjfcle,  since  it  yields 
a degenerate  distribution  at  aero  which  being  truncated  at  aero  does  not 
exist.  But  we  can  always  choose  A to  be  very  snail,  say  A * 0.00001, 
and  uae  this  value  of  A as  an  estimate  of  A.  This  approach  is  also 
supported  by  the  plot  of  the  norma  lined  likelihood  function  in  Figure 
3.2  for  this  deta  set. 

A similar  problem  may  arise,  when  we  try  to  find  the  moment 
estimators  or  the  modified  moment  estimators  proposed  by  Brass  (1958) , 
which  are  used  as  initial  estimates  far  the  maxiiiun  likelihood  estimation. 
Che  of  the  following  cases  may  arise: 

(i)  A and  y are  both  positive 

A A 

(ii)  A and  y are  both  negative 

(iii)  A or  y Is  equal  to  zero 

CSaea  (i),  (ii),  and  (iii)  will  occur  if  x < S2,  x > S2  and  x ■ S2, 
respectively.  If  x > S2  or  if  A and  y are  negative,  the  truncated  bi- 
nomial distribution  will  give  a better  fit.  But  here  we  are  considering 
the  distribution  of  the  lumber  of  failures  which  is  a truncated  ocnpound 
Poisson  distribution  and  has  a population  variance  _>  the  population  mean. 
In  such  a case,  the  truncated  Poisson  distribution  with  its  parameter  A 
obtained  from 

x - — will  be  a good  choioe. 
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Wok  the  sake  of  aomparison,  mb  also  obtained  estimates  of  y etc. 
using  Brass's  method  of  modified  moments.  The  results  for  the  eight  data 
sets  are  tabulated  in  Table  5.2.  A ocqpariaon  of  the  results  in  *»>hW 
5.1  and  5.2  sham  that  the  variability  of  Brass's  estimators  is  much 
larger  than  that  of  the  MLE's.  Also,  as  pointed  out  in  Section  4.2.3, 

the  A.R.E.  of  Brass's  estimators  for  the  range  of  parameters  of  interest 
here  is  about  75%,  a rather  poor  value. 

The  plots  of  the  fitted  inverted  gamma  densities  and  emulative 
distribution  functions  axe  given  in  Figures  4.1  to  4.8. 


Fig.  3.1.  Plots  of  the  Nbtrnalizad  Likelihood  Functions  For  the  Data  in  Table  1.2 

JM). 125  (.125)  2.0 


TABLE  4.1 


j^jfeaMINMTON  OF  PRIOR;  MTI  REFLECTOR  SM-225 

( Y-  679.0,  X-  0.2987,  Sf-  342.4,  o£-  0.3333,  P-  0.869) 
CtmrvA  and  Predicted  RmiuencleB 


1 

2 

3 

4 

5 

6 

7 

8 

9-19 


X Goodn— « of  fit  Tteet 


1 

2 

3-4 

5-7 

8-19 


n. 


13 

8 

8 

4 

8 


13 

8 

4 

4 

1 

2 

1 

4 

4 


41 


/V 

"x 


(Vfx) 


13 

7.3 

8.2 

5.9 

6.6 


0 

0.7 

-0.2 

-1.9 

1.4 


13 

7.3 

4.8 

3.4 

2.5 

1.9 

1.5 
1.2 
5.4 


41 


0 

0.067 

0.005 

0.612 

0.297 


0.981 


2 2 

x observed  “ °*98  » X 2,0.05  " 5,99>  81,3  P<X22  > X2obaervwd)  - 0. 
Thersfow,  truncated  negative  binomial  is  a good  fit. 
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TAHUE  4.2 


EETEWCNATION  or  PRIORI  SEARCH  MVPS  FP-3132 
( v-  713.7,  0.8646,  o*»  210.6,  0.3065,  P-  0.892) 


Obaarvad  and  Predicted  Frequencies 


A 

/\ 

X 

*Sc 

"x 

X 

"x 

nx 

1 

11 

12.4 

11 

0 

2.0 

2 

11 

10.0 

12 

0 

1.7 

3 

>10 

8.2 

13 

1 

1.5 

4 

7 

6.8 

14 

3 

1.2 

5 

5 

5.6 

15 

1 

1.1 

6 

4 

4.7 

16 

3 

0.9 

7 

1 

4.0 

17-44 

3 

5.3 

8 

4 

3.4 

9 

7 

2.8 

10 

3 

2.4 

””53 

>4 

74 

y2  Goodness  of  Fit  Ttest 

X 

"x 

Sc 

1 

ll 

12.4 

-1.4 

0.158 

2 

11 

10.0 

1.0 

0.100 

3 

10 

8.2 

1.8 

0.395 

4 

7 

6.8 

0.2 

0.006 

5 

5 

5.6 

-0.6 

0.064 

6-7 

5 

8.7 

-3.7 

1.574 

8-9 

11 

6.2 

4.8 

3.716 

10-12 

3 

6.1 

-3.1 

1.575 

13-44 

11 

10.0 

1.0 

0*  100 

• 

7.688 

2 

x observed 

7.688  and  »2  , „ M 

- 12.6, 

»<x2* 

Therefore,  tnwoatad  negative  binomial  is 

a good 

fit. 

) - 0.26 


obaarvsd  " 9,392  ****  4f0.0S  “ 9,49,»  “**  ^(xj  > X, 
*■**•»  truncated  nagativ*  binomial  Is  a good  fit. 
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TABUS  4.4 


CEfZEFMINATIQN  CF  PRIOR:  VIDEO  AMPLIFIER  AM-472 


( Y-  1581.5,  0.6675,  0^-779.7,  0.5277,  P*  0.923) 


A 

X 

*Sc 

1 

18 

17.7 

2 

11 

10.8 

3 . 

6 

7.0 

4 

6 

4.7 

5 

1 

3.2 

6 

2 

2.2 

7 

3 

1.6 

8 

1 

1.1 

9 

1 

0.8 

10-12 

2 

1.9 

51 

51 

x2  Goodneee  of  Tit  Test 

. 

X 

* 

, A . 2 .A 

W '"x 

1 

18 

17.7 

0.3 

0.005 

2 

U 

10.8 

0.2 

0.004 

3 

6 

7.0 

-1.0 

0.143 

4-5 

7 

7.9 

-0.9 

0.103 

6-12 

9 

7.6 

1.4 

0.258 

0.513 

*2  <tam  ’ °-5U  ■*  »2 

2,0.05  “ 

5.99,;  and  P(X22 

51  ^observed* 

Therefore,  truncated  negative  binomial  ia  a good  fit. 
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TABLE  4.5 


ULUMtCNATICN  OF  PRIOR: 


SYNCH  POWER  SUPPLY  SN-315 


( 7-  1246.0,  \m  0.00001,  - 947.4,  o%  - 0.4391,  p - 


Obearved  and  Predicted  Frequencies 


X 

"x 

A 

n 

X 

1 

22 

20.9 

2 

7 

8.0 

> 3 

1 

4.1 

4 

4 

2.4 

5 

2 

1.5 

6 

1 

1.0 

7 

3 

2.1 

40 

40 

X2  Goodness  of  Fit  Itest 

X nx 

A 

*x 

'vty 

1 22 

20.9 

1.1 

2 7 

8.0 

-1.0 

3-4  5 

6.5 

-1.5 

5-7  6 

4.6 

1.4 

TABLE  4.6 


CElEFMZNKnGN  OF  PRIOR:  SEARCH  INDICATOR  IP-128A 


( y-  1151.8,  A-  0.4407,  a~  542.0,  o*=  0.3993,  P»  0.905) 
Observed  and  Predicted  Frequencies 

x "x 


1 21 

2 8 

3 8 

4 2 

5 5 

6 4 

7 3 

8 2 

9-19  2 


19.4 

11.0 

7.1 

4.8 

3.4 

2.4 

1.8 
1.3 
3.8 


55 


55 


x2  Goodness  of  Fit  Teat 


X 

nx 

<V\> 

<VSt>2^ 

1 

21 

19.4 

1.6 

0.13 

2 

8 

11.0 

-3.0 

0.82 

3 

8 

7.1 

0.9 

0.11 

4-5 

7 

8.2 

-1.2 

0.17 

6-19 

11 

9.3 

1.7 

0.31 

1.54 


**  otwrvrt  - 1'54  and  *2  2,0.05  ’ 5-M ' p<X22  » *2olJ«rw^,  * 
Therefore,  truncated  negative  binomial  is  a good  fit. 


' tjOWS  4 ewjfta 
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TABI£  4.7 

ngflswnwnow  of  prior;  amft.tfier  panel  pn-8 

A ^ 

( Y-  754.9,  X-  0.4808,  o*.  309.9,  0.3286,  0.891) 

Observed  and  Predicted  Frequencies 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13-15 


Oooftiess  of  Fit  Tteat 


1 

2 

3 

4 

5-6 

7-9 

10-15 


58 

W 

15.8 

2.2 

10.0 

0.0 

7.0 

-2.0 

5.2 

-3.2 

7.1 

-1.1 

5.9 

1.1 

7.0 

3.0 

<v£>H 


0.306 

0.000 

0.571 

1.969 

0.170 

0.205 

1.286 

4.507 


* obeerved  " 4*507  «*  X2  4#0.05  - 9.49  , and  P(x24  > X2ob8erV6d) 
7ber»fore,  truncated  negative  binomial  is  a good  fit. 


TABUS  4.8 


flflliswfittanCN  OF  PRIOR:  TRftMSMTnry 


( Y*  680.3,  X»  0.2054,  <r^m  317.7,  oj«  0.2794,  p»  0.856) 
ObflBTwsd  and  Predicted 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10-33 


16 

10 

7 

5 

3 

3 

1 

1 

1 

3 


50 


x Goodneaa  of  Fit  tet 
x 


n 


1 

2 

3 

4-5 

6-8 

9-33 


16 

10 

7 

8 
5 
4 


nx 


17.5 

9.1 

5.8 

4.0 

2.9 

2.2 
1.7 
1.3 

1.0 

4.5 


50 


/V 

nx 

'v5*1 

17.5 

-1.5 

0.129 

9.1 

0.9 

0.089 

5.8 

1.2 

0.248 

6«9 

1.1 

0.175 

5.2 

-0.2 

0.008 

5.5 

-1.5 

0.409 

’ 1-°“  *2  3,0.05  - 7M<  *<*21  > X2, 

Ther*foz»,  truno«t«l  nagatiw  bimaial  la  a good  fit. 


1.055 

) - 0.79. 
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UtBUB  5.1 

PAWtSO^j*  ESTHETES  AW)  STWDMP  ERRORS  USING  MBBOP 


or  MMJm>  T.TK^r.TwnrQ 


Squipnant 

A 

Y 

Sy 

A 

A 

A 

P 

1. 

MCI  Fcf lector 
SH-225 

679.0 

342.4 

0.2987 

0.3333 

0.869 

2. 

Search  MVPS 
PP-3132 

713.7 

210.6 

0.8646 

0.3065 

0.892 

3. 

Oscilloscope 

06-126 

651.7 

256.1 

0.3335 

0.2616 

0.857 

4, 

Video  Amplifier 

MI-472 

1581.5 

779.7 

0.6675 

0,5277 

0.923 

5. 

Synch  Power  Sipply 
SN-315 

1246.0 

947,4 

0.00001 

0.4391 

0.908 

6. 

Search  Indicator 
ZP-128A 

1151.8 

542.0 

0,4407 

0.3993 

0.905 

7. 

Amplifier  Panel 
PH-8 

754,9 

309.9 

0.4808 

0.3286 

0.891 

8. 

Tranamitter 

680,3 

317.7 

0,2054 

0,2794 

0.856 
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TABUS  5.2 


PARAMETER  ESTIMATES  AND  STMCABD  ERRORS  USING 
BRASS'S  MODIFIED  METHOD  OF  MOMENTS 


Bguipaent 

A 

Y 

A 

a9 

A 

\ 

A 

A 

P 

1. 

MTT  Reflector 
34-225 

775.9 

2742.8 

0.3838 

2.5321 

.8943 

2. 

Search  MVPS 
pp-3132 

619.7 

1808.9 

0.7508 

2.6950 

.9033 

3. 

Oacilloeoope 

OS-126 

523.4 

1861.9 

0.2781 

2.0944 

.8764 

4. 

Video  Jtaplifier 
AM-472 

1754.0 

6529.5 

0.7696 

4.2833 

.9323 

5. 

synch  Power  Supply 
SN-315 

1676.1 

7832.9 

0.1290 

3.2113 

.9134 

6. 

Search  Indicator 
IP-128A 

1134.2 

4287.4 

0.3915 

2.9872 

.9074 

7. 

Amplifier  Panel 

1058.5 

3445.2 

0.7374 

3.3014 

.9161 

8. 

Tranenittar 

484.8 

1942.2 

0.0795 

1.8580 

.8611 
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far  Data  Set  2 


Figure  4.4  Plots  of  the  Fitted  Inverted  Gene  pdf  g(8)  md  odf  G(6) 
for  Data  Set  4 (y  - 1581.5,  A - 0.6675) 


Figure  4,6  Plots  of  the  Fitted  Inverted  Gamma  pdf  g(0)  and  cdf  G(8) 
for  Data  Set  6 (y  = 1151.8,  X - 0.4407) 
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6.  OnCXUPPC  ROBVRICS 


In  thin  report  m presented  a discussion  of  the  implications  of  the 
existence  of  a prior  distribution  for  the  parens  ter  of  the  exponential 
failure  distribution.  Tm  models  were  dismissal  to  describe  the 
heterogeneity  in  the  product. 

The  oaee  when  the  perneeter  6 has  an  inverted  ganma  density  was 
discuaasfl  in  detail.  The  method  of  obtaining  max  bran  lilcelihood  estimates 
for  the  parmasters  of  the  inverted  ganma  distribution  was  discussed  for 
tbe  cats  then  the  data  are  available  for  the  nmfcers  of  failure  in 
fixed  testing  times.  TVo  cases  mere  discussed:  the  ocnpleta  negative 
bincmial  and  the  truncated  negative  binomial. 

Eight  sets  of  field  failure  data  from  Schafer  (1970)  were  analysed 
in  detail  to  obtain  the  appropriate  fitted  distribution  and  the  results 
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APPENDIX  A 


SUBKVnW  1NBZM0 

PUBPOSE  AN?  PRXS3UHE: 

Ths  purpose  of  this  SUBROUTINE  is  to  fit  a truncated  negative  binomial 
distribution  to  the  observed  failure  data. 

For  the  given  data,  the  initial  parametric  estimates  are  obtained 
using  the  method  due  to  Brass.  Then  the  lmodiam  likelihood  estimtes  are 
obtained  based  on  the  expressions  given  in  Section  4.2  of  this  report. 

Also,  the  varianoe-aovarianoe  matrix  of  the  estimated  parameters  is 

calculated  using  the  expressions  of  Section  4.2. 

\ 

USAGE: 

CAUL  TNBINO  (N,T,F) 

PARAMETERS 

N - maxisun  nunbsr  of  failuree  in  the  data  set 
T - testing  time 

F - input  vector  of  failure  frequencies  (vector  If  length  N) 
OUTPUT: 

Initial  estimates  of  Goema  and  Lamia. 

Variance  coverisnos  matrix  of  - tj—<» 

Correlation  coefficient  between  rv— and  nmmm 

Table  of  obeerved  frequencies,  expected  frequencies  baaed  on  initial 

estimates  and  aaqpcted  frequencies  baaed  on  final  estimates. 

LISTING: 

A listing  of  the  SONOUnMB  is  given  on  the  following  page. 


SUBROUTINE  TNBINQ(N»T »F) 

TRUNCATED  NEGATIVE  BINOMIAL 
DIMENSION  X<50>  »F<50>»FI<50> ,AX<50>  »FFI<50> 

DIMENSION  YYY < 500 ) » ZZZ < 500) 

DO  900  lUU^lflO 
NO  I - 1 

DO  500  LBJ  = lrNOI 
PRINT  200 
200  FORMAT (1H1) 

300  FORMAT</> 

SF  - 0. 

SFX  = 0. 

SQ  * 0. 

DO  350  I « 1»N 
350  X(I)  - I 

C****COMPUTE  XBAR  AND  SQR 
DO  2 I = IfN 

2 SF  = SF+F< I ) 

DO  3 I = 1»N 

3 SFX  » SFX+X ( I ) #F  < I > 

XBAR  ■ SFX/SF 

DO  4 I * 1»N 

4 SQ  = SQ+F (I)#(X(I) -XBAR ) ##2 
SQR  ® SQ/(SF-1 .0) 

PRINT  101 » XBAR r SQR 

101  FORMAT <2X»6HXBAR  =»F10.4r5Xr4HS**2  =>F10.4/> 

C****COMPUTE  INITIAL  ESTIMATES  BY  METHOD  DUE  TO  BRASS 
Z - XBAR#T# ( SF-F ( 1 > ) / < SF#SQR~XBAR# ( SF-F ( 1 ) ) ) 

ZZ  * < XBAR* (Z/<T+Z)>-(F<1) /SF )>/(l«~<Z/( T+Z ) ) > 

PRINT  25rZrZZ 

25  FORMAT <2X»7HGAMMA  = » F20 ♦ 5 » 5X  r 8HLAMBDA  =»F10.5r/> 

FI  ( 1 ) = (ZZ*(<Z/n+Z))**ZZ)*(T/<T+Z)))/<l.-<<Z/(T+Z))4t*ZZ)) 
FI < 1 ) » SF*FI < 1 ) 

DO  20  I = 2»N 

20  FI < I ) ■ <<ZZ+I- 1)/I)*<T/<T+Z))*FI<I-1> 

C****MAXIMUM  LIKELIHOOD  ESTIMATES 
Y ■ Z 
Z ■ ZZ 
AX< 1 ) ■ SF 
DO  5 I ■ 2*N 

5 AX  ( I > » AX(I-1)--F(1-1) 

DO  400  JFK  » 1»20 
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DO  400  JFK  » If 20 
S = 0. 

SS  a o. 

DO  6 I = 1 f N 

S = S+<1./  < Z+I-l ) ) *AX < I ) 

A SS  a SSMl./<  <Z+I-1  >*#2>)*AX<I) 

£i  a Y/<T+Y) 

£ a l.-D 

DZ  » D**Z 

A a SF#Z#(E/Y)/< 1 »-DZ)-8FX/<T+Y) 

B = SF'*Al.OG  < D ) / < 1 ♦ ~DZ ) +S 

CC  =-  SF* < E/Y ) ♦ < 1 » -DZ# < 1 ♦ -Z#ALOG (D)))/<(1* -DZ > ##2 ) 

AA  = SFX/<  <T+Y>#*2)  - <SF#Z*T/<Y#Y#<T+Y)*<T+Y>  >)*<T-F2,*Y-DZ* 
1<T+2.*Y+Z*T>)/<<1.-DZ>##2> 

BB  * SF#DZ*<  (ALOG(D) )##2>/< ( 1 ♦-DZ)##2)-SS 

GG  a AA*BB-CC*CC 

PRINT  7fZf AfBf AArBBfCCfGG 

7 FORMAT ( 2X  f 8E1 2 ♦ 5 ) 

YD  = - < BB/GG ) *A+  < CC/GG ) #B 
ZD  ==  ( CC/GG  )*A~<AA/GG)#B 
Y a y+YD 
Z a z+ZD 

IF ( ABS < A ) -0 ♦ 001 ) 8f 8f 400 

8 1F(ABS(B)-0*001 ) 9f9f400 
400  CONTINUE 

C#*#*OARIANCE~COOARIANCE  MATRIX  OF  MLE 

9 Y=Y-YD 

Z a Z-ZD 
PRINT  300 
PRINT  25r YfZ 
YYY(LBJ)  a y 
ZZ7.(LBZ ) « Z 
OG  « -BB/GG 
SOG  a SORT < 00) 

PR  IN!  12 f OG f BOG 

12  FORMA T < 2Xf  19M0ARIANCE  OF  GAMMA  * f F20 . Sif  19HSTD . DEO . OF  GAMMA  a, 
ir20«[if/) 

01.  ; -AA/UG 

sol  a sn«r(OL> 

i f I NT  13f0LfS0L 

13  F ORMAT < 2X r 21HVARI ANCE  OF  LAMBDA  » rF20.3»20HSTB.DEU.  OF  LAMBDA 
ll-'20.5f/) 

OLG  = CC/GG 
PRINT  14  f OuG 

14  F ORMAT <2Xf25HCD0AR I ANCE  LAMBDA- GAMMA  arF20»5r/) 
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14  FORMAT <2Xr25HCQUARIANCE  LAMBDA-GAMMA  *rF20.5»/> 

RHO  - ULG/SQRT  < VL4VG  > 

PRINT  15f RHO 

15  FORMAT <2X»5HRH0  =»F20.5r/> 

C******CALCULATIONS  IN  TERMS  OF  MEAN 

XM  - Z*T/Y 
PRINT  30rXM 

30  FORMAT <2X»6HMEAN  =rF20.5f/) 

XA  = Z/CXM+Z) 

AAA  = XA*<-SFX/(XM*XM>+SF*<XA**<Z+1.))/<(1.-XA**Z>#*2>) 

CCC  '»  -SF*<XA**<Z+1,)>*<1.-XA+AL0G<XA>>/<<1,-XA**Z>**2> 

BBB  - BB+CCC*CCC/AAA-CC*CC/AA 
GGG  * AAA*BBB-CCC*CCC 
UM  ■ -BBB/GGG 
SVM  * SORT  < VM) 

PRINT  31 » VM»SUM 

31  FORMAT <2X»18HUARIANCE  OF  MEAN  «,F20.3r 18HSTD.DEV.  OF  MEAN  =*F20. 

It/'t 

WL  -AAA/GGG 
SVVL  - SORT ( WL) 

PRINT  13»WL»SWL 
VLM  * CCC/GGG 
PRINT  32rVLM 

32  FORMAT (2Xr24HC0UARIANCE  LAMBDA-MEAN  =,F20.5,/> 

RRHO  * ULM/SCJRT  < UVL4VM  ) 

* PRINT  15 » RRHO 

C*****************t******t****^***********M*****t******t*****t^****** 

FFI  < 1 ) - SF*Z*DZ*E/<i.~DZ) 

DO  16  1 » 2 r N 

16  FFI < I ) * (<Z+I-1)/I )#E*FFI ( 1-1 ) 

17  FORMAT <5X » 1HX  r 12X * 1HF  1 14X»2HFI » 14X  >2HFI ) 

DO  18  I - lrN 

PRINT  19»X(I)»F<I)»FI(I)»FFI(I) 

19  FORMAT ( 3X » F4 • 0 » 8X » F5 ♦ 0 r 5X  * F 1 2 ♦ 4 » 5X » F 1 2 ♦ 4 > 

18  CONTINUE 
500  CONTINUE 
900  CONTINUE 

RETURN 

END 


*1 


